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Incorporation  Kinetics  in  Mixed  Anion  Compound  Semiconductor  Alloys 

Joanna  M.  Millunchick  and  Evan  M.  Anderson,  Department  of  Materials  Science  and 
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Abstract 

We  present  a  kinetic  modei  predicting  anion  incorporation  in  inAsSb.  inciuded  are 
the  effects  of  As  desorption,  Sb  segregation,  and  Sb  dispiacement  by  As,  any  of  which 
may  be  iimited  by  the  in  flux  if  it  is  comparativeiy  iarger.  The  modei  captures 
experimentai  data  over  a  range  of  growth  conditions  for  the  inAsSb  system  using 
activation  energies  for  desorption  and  Sb  segregation  found  in  iiterature.  The 
activation  energy  for  Sb  dispiacement  found  in  this  work  is  l.SeV.  This  modei  is 
generai  and  shouid  be  vaiid  for  other  mixed  anion  systems,  or,  appropriate^ 
modified,  mixed  cation  systems  and  mixed  anion/cation  systems  such  as  AiinAsSb. 
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Since  the  invention  of  the  MBE  technology,  most  work  on  alloys  with  three  or 
more  constituents  has  focused  on  mixed  Group-111  materials.  This  has  enabled  the 
development  of  a  wide  range  of  important  semiconductor  heterostructures  and 
devices.  Alloying  the  Group-111  elements  was  a  natural  beginning  for  the  technology 
due  to  the  relatively  straight-forward  composition  control  allowed:  under  most 
typical  growth  conditions  the  Group-111  elements  add  linearly  to  the  total  growth 
rate  and  the  resulting  composition  is  directly  obtained  as  a  ratio  of  the  individual 
rates.  In  contrast,  mixed  anion  (Group-V)  alloys  are  more  complex  to  control.  The 
Group-V  elements  usually  do  not  influence  the  growth  rates,  and  the  resulting 
composition  is  not  merely  the  ratio  of  the  fluxes.  Instead,  there  are  many  competing 
processes  that  alter  the  incorporation  rates  of  these  elements.  However,  mixed 
Group-V  alloys  offer  many  new  important  materials  properties  that  are  not  possible 
with  Group-111  alloys  and  are  therefore  of  interest  to  explore.  Examples  are  dilute 
bismuthides  ^  and  nitrides  ^  that  can  produce  reduced  bandgaps  while  remaining 
lattice  matched  to  available  substrates,  and  the  emerging  field  of  highly-mismatched 
alloys  such  as  dilute-antimonide  GaNSb  ^  and  very  narrow  bandgap  materials  such 
as  InAsSb,  which  is  the  main  focus  of  this  work. 

Mixed  As/Sb  alloys  have  been  produced  for  several  decades.^  The  properties 
of  MBE-grown  InAsSb  in  particular  has  been  studied.  but  few  reports  attempting 
to  systematically  determine  the  Sb/As  incorporation  as  a  function  of  growth 
parameters  exist  other  than  some  empirical  work.  It  has  recently  been 
demonstrated  that  InAsSb  has  a  smaller  bandgap  than  previously  thought.  This 
opens  up  a  long  sought  path  to  a  111-V,  direct-bandgap,  long-wavelength,  infrared 
detector  material  that  can  rival  HgCdTe  without  the  complications  and  drawbacks 
associated  with  superlattices. In  order  to  produce  such  material,  it  needs  to  be 
grown  strain-free  so  that  the  resulting  bandgap  is  not  affected  by  residual  strain  or 
Group-V  sublattice  ordering,  and  the  composition  of  the  alloy  needs  to  be  well 
controlled  so  that  the  bandgap  targeted  for  specific  device  applications  is 
obtained. 14  Further  materials  optimization  requires  variations  of  the  growth 
conditions,  while  maintaining  a  desired  target  composition.  It  is  therefore  of  great 
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importance  to  be  able  to  model  the  incorporation  process  to  guide  and  reduce  the 
number  of  optimization  experiments. 

Typically,  any  general  kinetic  model  takes  into  account  all  of  the  physical 
processes  that  take  place  during  growth.  Many  such  models  have  been  developed 
for  the  mixed  anion  alloys,  including  GaAsP  GaAsSb  and  GaAsBi  No 
modeling  of  the  growth  of  InAsSb  is  available  in  the  literature,  however.  Here  we 
present  a  model  that  addresses  incorporation  in  unstrained  material  taking  into 
account  desorption,  Sb  segregation  and  As  displacement  of  Sb  and  apply  it  to  InAsSb. 
The  only  fitting  parameter  of  the  resulting  model  is  the  activation  energy  of  As 
displacement,  which  was  found  to  be  1.3  eV.  The  model  fits  agree  well  with 
experimental  results  obtained  in  multiple  laboratories,  and  provide  the  sought  after 
ability  to  predict  the  Sb  composition  from  known  experimental  parameters. 

Kinetic  model 

In  this  general  kinetic  model,  we  consider  the  incorporation  of  elements  into 
a  compound  MAi-xSx,  where  M  is  the  Group-111  metal,  A  is  the  Group-V  anion,  and  S  is 
the  Group-V  segregant,  that  is,  the  element  that  is  more  easily  displaced  in  the  alloy 
(eg:  Sb  in  InAsi-xSbx  or  Bi  in  GaAsi-xBix).  We  are  interested  in  understanding  how 
the  concentration  of  the  segregant  in  the  bulk  Xb  changes  as  a  function  of  growth 
parameters  such  as  temperature  T,  flux  of  the  segregant  Fs,  flux  of  the  anion  Fa,  and 
flux  of  the  metal  Fm.  We  do  this  by  taking  into  account  all  of  the  processes  that  affect 
the  coverage  of  each  of  the  constituents.  The  composition  of  the  segregant  at  the 
surface  Xs  is  a  ratio  of  the  coverage  of  each  species  in  the  following  manner: 


Xs  =  0s  /(0s  +  0a) 


(1) 


where  0  is  the  surface  coverage  of  the  anion  and  segregant  respectively.  Xb  and  Xs 
are  not  necessarily  the  same  value  due  to  surface  segregation,  as  will  be  discussed 
below. 
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Several  processes  have  been  identified  that  affect  the  amount  of  each  species, 
as  shown  in  Fig.  1.  The  incident  flux  is  the  primary  term,  but  this  contribution  is 
reduced  in  the  presence  of  desorption.  Surface  segregation,  or  the  exchange  of  a 
surface  anion  species  with  a  subsurface  segregant  species,  is  also  well  known  for  a 
number  of  mixed  anion  alloys  such  as  GaAsSb,  GaSbBi,  and  inAsSb.  Both  Sb  and  Bi 
have  been  reported  to  be  excellent  surfactants,  remaining  primarily  on  the  surface 
rather  than  incorporating  into  the  film  under  the  appropriate  growth  conditions. 
i^it  has  also  been  reported  in  some  mixed  anion  alloys  that  one  species  may  displace 
the  other.  For  instance,  it  has  been  shown  that  the  incorporation  of  Sb  in  GaAsSb  is 
hampered  in  the  presence  of  a  large  As  flux,  ^  however,  this  effect  has  not  been 
reported  for  inAsSb. 

The  rate  of  desorption,  displacement,  and  segregation  are  activated 
processes,  and  have  a  strong  dependence  on  temperature  with  the  standard 
Arrhenius  form: 

Pj  =  iVvexp(-^)  (2) 

Where  N  is  the  number  of  atomic  sites  on  the  surface,  v  is  the  attempt  frequency 
4x10^2  s-i  ^  gj  is  i;he  activation  energy  of  the  ith  process,  k  is  Boltzmann’s  constant, 
and  T  is  the  temperature,  it  is  also  possible  to  include  terms  for  the  formation  of 
metal  or  segregant  droplets,  which  has  been  observed  in  growth  with  Sb  and  Bf^i 
but  for  the  purpose  of  this  work  we  will  only  consider  anion  overpressures  for 
which  droplets  are  not  expected  to  form,  in  this  case,  the  surface  is  always  covered 
in  either  the  anion  or  segregant  such  that  0s  +  0a  =1. 

The  rate  of  each  process  is  not  the  only  factor  that  needs  to  be  included.  The 
form  of  each  term  must  also  include  a  limiting  factor  (t)ij  that  takes  into  account  the 
fact  that  deposition  of  another  species,  Fj,  may  limit  the  process  in  question: 

^0-=?^  (3) 

(t)ij  approaches  unity  when  the  reaction  rate  is  fast  compared  to  the  arrival  of  the 
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species  that  would  prevent  the  reaction  from  taking  place,  and  zero  when  the 
reaction  rate  is  slow  in  comparison.  In  the  experiments  that  follow,  the  limiting  flux 
is  always  that  of  the  metal,  thus  we  simplify  the  notation  for  the  limiting  factor  to  (j)!. 


The  fact  that  there  is  surface  segregation  implies  that  there  is  a  difference 
between  the  composition  at  the  surface  Xs  and  in  the  bulk  Xb.  The  bulk  composition 
is  proportional  to  the  surface  composition  via  the  equilibrium  partition  coefficient  K, 
which  is  the  fraction  of  the  segregant  that  remains  in  the  film. 22  The  value  of  K  can 
range  from  0  (for  a  pure  surfactant)  to  1  (for  complete  incorporation).  That 
coefficient  may  be  modified  via  the  ratio  of  the  rate  of  the  segregation  Pseg  and  the 
rate  of  the  arrival  of  the  metal,  Fm,  that  frustrates  the  segregation.  The  form  is  such 
that  the  effective  coefficient  goes  to  K  when  Fm  is  small  compared  to  Pseg,  and  unity 
when  it  is  large.  Thus  the  composition  of  the  bulk  is  the  product  of  the  surface 
concentration  and  the  effective  segregation  parameter: 


(4) 


The  growth  rate  of  the  segregant  is  given  by: 

Js  ~  Fs  ~  6sPsdes4*sdes  ^b6aPseg4*seg  ~  6sPdis4*dis  (S) 

The  first  term  on  the  right  hand  side  represents  the  amount  of  impinging  material, 
Fs.  The  second  term  represents  the  reduction  in  the  net  flux  due  to  desorption, 
which  is  the  product  of  the  rate  of  desorption  and  the  limiting  factor  comparing  the 
desorption  rate  and  the  arrival  of  the  metal.  The  third  term  increases  the  net  flux  via 
surface  segregation,  which  is  the  product  of  the  rate  of  segregation  and  the  limiting 
factor  comparing  the  desorption  rate  and  the  arrival  of  the  metal.  The  fourth  term 
reduces  the  net  flux  as  a  result  of  anion  displacement,  which  is  the  product  of  the 
rate  of  displacement  and  the  limiting  factor  comparing  the  desorption  rate  and  the 
arrival  of  the  metal.  The  growth  rate  of  the  anion  species  may  be  written  in  a  similar 
manner,  with  the  exception  of  terms  3  and  4,  which  describe  reaction  rates  between 
the  anion  and  segregant.  The  displacement  term  is  the  same  for  both  equations. 
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except  that  this  term  acts  to  increase  the  anion  growth  rate.  The  segregation  term, 
on  the  other  hand,  acts  to  decrease  the  anion  growth  rate.  These  equations  can  be 
solved  simultaneously  to  determine  Xb  as  a  function  of  the  growth  parameters. 

Experimental  Data 

Several  InAsi-xSbx  samples  were  grown  with  two  MBE  systems  in  our 
respective  laboratories  and  compared  to  the  predictions  of  the  kinetic  model.  All 
samples  were  growth  using  solid  sources  for  Group-111  and  valved  sources  with 
cracking  zones  for  As  and  Sb.  ARL  samples  were  grown  with  a  residual  strain  less 
than  0.1%  on  graded  buffer  layers  as  described  earlier  ^3  at  a  substrate  temperature 
of  415  C,  measured  with  a  K-space  BandlT  system  operated  in  pyrometry  mode  and 
calibrated  with  a  blackbody  source  through  the  heated  viewport.  The  emissivity  was 
set  to  produce  the  characteristic  RHEED  pattern  transition  on  GaSb  from  x3  to  x5  at 
430  C.  24 

Absolute  flux  values  were  determined  in  two  different  ways.  In  the  first  case 
we  started  by  measuring  the  growth  rates  with  RHEED  oscillations  under  excess 
group  V  flux,  which  allowed  us  to  relate  the  group  111  BEP  to  flux.  We  then  reduced 
the  group  V  flux  until  we  observed  a  reduction  in  the  growth  rate,  at  which  point  we 
assumed  that  the  group  111  and  V  fluxes  were  equal.  In  the  second  case,  flux  was 
determined  by  monitoring  the  RHEED  pattern  until  it  changed  from  a  group  V 
terminated  surface  to  a  group  111  terminated  surface.  This  point  was  taken  to 
correspond  to  a  V/111  ratio  of  1.  Because  the  group  111  growth  rate  is  not  influenced 
by  desorption  or  any  other  depleting  processes,  the  RHEED  calibrated  rate  can  be 
assumed  to  be  equivalent  to  its  flux  as  measured  by  a  Beam  Flux  Monitor.  After 
taking  into  account  beam  flux  monitor  sensitivities,  the  In  flux  is  assumed  to  be 
equal  to  the  anion  flux  at  V/lll=l.  The  procedure  was  done  for  both  anions.  UM 
samples  were  grown  over  a  range  of  In  rates  using  different  temperatures,  as 
measured  using  a  low  temperature  pyrometer.  Film  growth  was  monitored  in  situ 
using  RHEED  and  ex  situ  using  high-resolution  x-ray  diffraction.  The  Sb- 
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composition  values  were  obtained  from  XRD  data  using  symmetric  (004)  scans  and 
asymmetric  (115)  scans. 

Figure  2  shows  the  bulk  composition  Xb  of  nearly  lattice  matched 
InAsSb/GaSb  as  a  function  of  growth  temperature  T  and  In  growth  rate  Fm.  The  data 
exhibits  a  small  negative  slope  for  each  set  of  data  at  a  constant  growth  rate,  and 
increasing  bulk  composition  with  increasing  Fm.  Figure  2  also  shows  the  fit  to  the 
data  using  the  kinetic  model  using  activation  energies  for  the  desorption  and 
segregation  found  in  the  literature  (shown  in  Table  1).  The  segregation  energy 
calculated  from  reference  is  1.2  eV,  and  the  segregation  parameter  is  0.4.  The 
activation  energy  for  As  displacement  of  Sb  was  a  fitting  parameter,  and  was  found 
to  be  1.3  eV.  The  fit  shows  that  the  temperature  dependence  of  the  composition  is 
not  a  simple  Arrhenius,  but  instead  depends  on  the  interplay  between  the  different 
terms.  In  these  experiments,  the  desorption  rates  are  slow  and  don't  contribute 
much  to  the  fluxes  and  thus  to  the  composition.  Both  the  segregation  and 
displacement  terms  act  to  reduce  the  composition  as  long  as  Fm  is  low.  Since  the 
activation  energy  to  segregate  is  lower  than  to  displace  Sb,  segregation  is  the 
dominant  mechanism  that  reduces  composition  at  low  temperature.  Assuming  only 
segregation  the  composition  will  increase  with  increasing  temperature  until  it 
reaches  the  equilibrium  composition,  0.4.  The  fact  that  lower  compositions  are 
observed  suggest  that  the  rate  of  displacement  of  Sb  by  As  becomes  appreciable  at 
higher  temperatures,  which  reduces  the  coverage  of  Sb  on  the  surface  and  thus  the 
amount  of  Sb  that  may  incorporate  into  the  film. 

The  fact  that  the  composition  increases  with  increasing  In  flux  is  evident 
from  the  data.  The  model  also  predicts  this  trend.  An  increase  in  Fm  increases  the 
amount  of  Sb  incorporation.  Inspection  of  equation  3  shows  that  an  increase  in  Fm 
will  result  in  a  decrease  of  the  efficiency  factor  (j)  for  both  the  As  displacement  and 
Sb  surface  segregation.  Both  of  these  effects  lead  to  more  Sb  incorporation. 
Furthermore,  higher  Fm  will  result  in  an  increase  in  the  segregation  coefficient  (equ. 
4).  All  of  these  indicate  that  increasing  Fm  will  result  in  the  kinetic  trapping  of  the 
segregant.  The  desorption  rates  for  As  and  Sb  are  slow  for  InAsSb,  ^6,27  typically 
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less  than  the  arrival  rate  of  In,  so  that  these  terms  have  a  negligible  impact  on  the  Fm 
dependence  of  the  composition. 

The  plots  in  Fig.  3  demonstrate  the  dependence  of  the  bulk  composition  as  a 
function  of  group  ill  and  V  fluxes.  Fig.  3a  shows  the  Sb/As  flux  ratio  at  the  growth 
temperature  of  415  C.  The  points  of  the  same  color  are  those  for  which  the  As  flux 
was  varied  and  the  in  and  Sb  were  held  nominally  constant.  The  experimental  data 
shows  that  the  composition  increases  sublinearly  with  increasing  Sb/As.  The  solid 
lines  show  the  fits  of  the  model  to  the  data  using  the  parameters  listed  in  Table  1. 
The  model  predicts  that  there  is  a  dependence  on  the  absolute  in,  Sb,  or  As  flux  on 
the  composition  as  a  function  of  Sb/As  flux  ratio,  in  agreement  with  results  shown 
for  GaAsSb.^  The  amount  of  Sb  incorporation  increases  with  lower  absolute  Sb  flux, 
presumably  because  the  As  displacement  is  lower.  The  model  also  predicts  that  the 
absolute  fluxes  influence  the  composition  as  a  function  of  V/iii  ratio,  shown  in  Fig. 
3b.  The  composition  decreases  with  V/iii  ratio  when  holding  the  Sb  and  in  fluxes 
constant.  The  magnitude  of  the  composition  depends  on  the  absolute  value  of  the  Sb 
(As)  flux,  with  higher  Sb  (As)  flux  resulting  in  a  higher  (lower)  incorporation  of  Sb. 
This  also  supports  the  observation  that  the  As  displacement  mechanism  is  offset 
(enhanced)  by  the  higher  flux  of  Sb  (As). 

Figure  4  shows  a  map  that  plots  the  growth  conditions  to  obtain  particular 
compositions  for  an  in  flux  of  6.3x10^^  corresponding  to  1  um/hr,  and  growth 
temperatures  of  405,  415,  and  425.  The  limit  on  the  y-axis  does  not  go  to  zero, 
because  ratios  less  than  1  indicate  that  the  surface  is  metal  terminated,  conditions 
for  which  this  model  is  not  valid.  The  contours  from  left  to  right  represent  the 
conditions  for  0.1  <x  <0.4.  The  composition  varies  weakly  as  a  function  of  the  V/iii 
ratio  at  low  Sb/As  ratio,  but  increases  in  sensitivity  as  the  Sb/As  ratio  increases. 
Similarly,  the  composition  is  weakly  affected  at  low  Sb/As  ratios  and  more  strongly 
at  higher  ratios.  Both  of  these  facts  arise  because  the  interaction  between  the  Sb 
segregation  and  the  As  displacement  becomes  stronger  as  the  absolute  amount  of 
both  anion  fluxes  increases. 
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Conclusions 

We  have  shown  that  a  relatively  simple  model  can  well  describe  the 
qualitative  and  quantitative  incorporation  of  As  and  Sb  in  InAsSb.  Both  As 
displacement  and  Sb  segregation  result  in  a  depletion  of  Sb  in  the  film.  The 
activation  energy  for  As  displacement  of  Sb  is  found  to  be  1.3  eV.  Any  of  these 
processes  may  be  limited  by  the  in  flux,  which  will  inhibit  that  process  if  it  is  larger 
than  the  process  rate,  however,  this  has  a  relatively  small  affect.  A  growth  condition 
map  was  also  generated  that  provides  valuable  guidance  on  how  to  achieve 
compositional  control  of  this  difficult  to  grow  alloy.  This  model  is  completely 
general  and  should  also  be  valid  for  other  mixed  anion  systems,  including  those  with 
phosphorous  or  bismuth.  Furthermore,  the  model  may  be  easily  modified  to 
examine  mixed  cation  alloys  where  one  cation  acts  as  a  segregant,  or  even  mixed 
anion/mixed  cation  systems  such  as  AlinAsSb. 
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Table  1:  Parameters  used  to  fit  the  model  to  the  experimental  data. 


Desorption  Energy  of  As  from  InAs 

2.4  eV 

Reference 

Desorption  Energy  of  Sb  from  InSb 

2.7  eV 

Reference  26 

Segregation  Energy  of  Sb  in  InAsSb 

1.2  eV 

Reference  26 

Equilibrium  segregation  parameter 

0.4 

Reference  25 

Energy  for  As  displacement  of  Sb  in  InAsSb 

1.3  eV 

this  work 
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segregaiit  aiiion  metal 


Figure  1:  Factors  governing  the  surface  coverage  include  (a)  deposition,  (b) 
desorption,  (c)  displacement,  and  (d)  segregation. 


Figure  2:  Sb  bulk  concentration  as  a  function  of  inverse  temperature  for  In  fluxes 
l<Fm<5.5xl0i^  cm-^s'i.  The  solid  lines  are  the  model  fit  using  the  parameters  in 
Table  1 
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V/lll  Ratio 


Figure  3:  Sb  concentration  as  a  function  of  (a)  Sb/As  ratio  and  (b)  V/III  ratio  at 
T=415  C  for  Fm=  6.3x10^18  cm-^s'i  (red  and  blue)  and  1.0x10^18  cm  ^s'i  (green). 
The  solid  lines  are  fits  to  the  model  using  the  parameters  in  Table  1. 
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Sb/As  Ratio 

Figure  4:  As  and  Sb  fluxes  required  to  obtain  certain  compositions  for  Fm=6.3xl0i^ 
cm-2s-i,T=405C  (dashed),  415C  (solid),  and  425C  (dotted). 
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